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ABSTRACT 

Context. We present a search for outflows towards 51 submillimetre cores in the Perseus molecular cloud. 

Aims. Our first objective is to identify the protostellar population through the detection of molecular outflows. Our second aim is to 
consistently derive outflow properties from a large homogeneous dataset within one molecular cloud in order to investigate further 
the mass dependence and time evolution of protostellar mass loss. 

Methods. We used the James Clerk Maxwell Telescope to map 2' x 2' regions around each core in '^CO 3-2. Where molecular 
outflows were detected we derived momentum fluxes. 

Results. Of the 51 cores, 37 show broad linewings indicative of molecular outflows. In 13 cases, the linewings could be due to 
confusion with neighbouring flows but 9 of those sources also have near-infrared detections confirming their protostellar nature. The 
total fraction of protostars in our sample is 65%. All but four outflow detections are confirmed as protostellar by Spitzer IR detections 
and only one Spitzer source has no outflow, showing that outflow maps at this sensitivity are equally good at identifying protostars 
as Spitzer. Outflow momentum flux correlates both with source luminosity and with core mass but there is considerable scatter even 
within this one cloud despite the homogeneous dataset. We fail to confirm the result of Bontemps et al. (1996) that Class I sources 
show lower momentum fluxes on average than Class sources, with a KS test showing a significant probability that the momentum 
fluxes for both Class Os and Class Is are drawn from the same distribution. 

Conclusions. We find that outflow power may not show a simple decline between the Class to Class I stages. Our sample includes 
low momentum flux, low-luminosity Class sources, possibly at a very early evolutionary stage. If the only mass loss from the core 
were due to outflows, cores would last for lO'-lO* years, longer than current estimates of 1.5-4 x 10^ years for the mean lifetime for 
the embedded phase. Additional mechanisms for removing mass from protostellar cores may be necessary. 

Key words. Submillimeter;Stars: formation;Stars: evolution;ISM: jets and outflows 



1. Introduction 

Submillimetre (submm) continuum surveys o f the Perseus 
molec ular cloud, such as our SCUBA survey , dHatchell et alJ 
I2005L hereafter Paper I, see also lKirk et al.l20d6l). a nd the similar 
1300//m survey with Bolocam ( Enoch et al.l 2006h . have identi- 
fied over 100 submm cores. These cores could be either proto- 
stellar or starless in nature: protostellar cores can be identified by 
radio or infrared detection of a centr al source, or evidenc e for a 
protostellar outflow. In a sister paper dHatchell et alj20d6l here- 
after Paper II) we classify the protostellar cores on the basis of 
their sp ectral energy distributions, including the new data from 
Spitzer dEvans et alj|2003l: |j0rgensen et al.ll2007h . Here we take 
a complementary approach and look for the presence of molec- 
ular outflows in order to separate the protostellar cores from the 
starless population. As outflow activity is thought to begin al- 
most as soon as the formation of a hydrostatic core, outflows are 
a good indicator of the presence of embedded protostars. 

In addition, studying molecular outflows is an important 
way to get information on the accretion/mass ejection process 
which drives protostellar evolution. Several important studies 
have shown that the outflow force or momentum flux Fco is 
corre lated with bo th mass and luminosity of the protostellar 
core dCabrit & Bert out 1992; Bontemps et al. 1996i) . a correla- 
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tion which r ecently has been shown to hold over several orders of 
magnitude dChurchweIllll997l: iBeutheret aP 120021: IZhang et all 
I2005h and held to be strong evidence for a common protostellar 
formation mechanism for high- and low-mass sources alike. The 
details of mass accretion and ejection, however, remain uncertain 
with competing mechanisms for the interaction between mass 
flow and magne ti c fields in the inner part s of accretion disks 
dShu et all [1991 IWardle & Koeiiigil [19931: iPelletier & Pu"drit3 
|1992 ) and considerable uncertainty on how the accretion rate 
varies with time. One certainty, however, seems to be that 
the accretion rate a nd the momentum output in the jet/wind 
are strongly linked (|Shu et al. I I 19941: iFerreira & Pelletieill 19951; 
iBontemps et al.lll996l) . and therefore, assuming a roughly con- 
stant mass entrainment fraction, one can use the outflow momen- 
tum as a surrogate measure of the momentum in the jet/wind, and 
therefore the mass accretion rate onto the protostar 

We therefore set out to make small outflow maps around the 
submm cores in Perseus to determine their outflow properties. 
With a consistent dataset for a moderately large sample drawn 
from one molecular cloud, we can look statistically at the out- 
flow detection rate and the variation of the momentum flux with 
other parameters. 

We assume a distance of 320 pc for the Pe rseus molecular 
cloud based on the Hipparcos distance of IC348 ( de Zeeuw et al.l 
1 19991) . consistent with our earlier paper (iPaper u) . though read- 
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ers should be awa re that other recent studies dKirk et al.l 120061 : 
lEnoch et"ari |2006|) have assumed a closer distance of 250 pc 
based on extinction studies ( ICernis & Straizvs 2003i) and refer- 
ences therein). 

In Sect. |2] we describe the observations. Sect.[3]contains the 
results and an examination of the detection statistics. Momentum 
fluxes and outflow evolution are considered in Sect.|4]. In Sect. [5] 
we discuss the derived outflow mass loss rates and the impli- 
cations for core destruction tim escales given the core masses 
implied by the submm emission (iPaper It iPaper It) . Finally, our 
conclusions are given in Sect.|6] 



2. '^CO 3-2 outflow observations 

Outflow maps were made in '^CO J = 3-2 in the years 
2000-2003 using RxB3 on the James Clerk Maxwell Telescope 
(JCMT). '^CO 3-2 is a better tracer of molecular outflows 
than lower-J transitions as it is more sensitive to warm outflow 
than cold ambient emission. The sample consisted of contigu- 
ous areas a round the known cluste rs Bl, IC348, L1448, L1455, 
NGC1333 dKnee & Sandelll |2000|) . and a number of isolated 
submm cores. The map areas were extended 2' x 2' around each 
isolated source or tiled to cover the cluster regions. The submm 
cores which lie within the outflow maps are listed in Table [1] 
P ositions and further i nformation on these cores can be found 
in lPaper j and lPaperlj The beamsize at 345 GHz is 14" and the 
maps were fully sampled at 5" spacing in RA and Dec. A correc- 
tion for main beam efficiency of t^mb - 0.63 was applied. The 
noise level was 0.5 K or less on 625 kHz (0.54 km s"') chan- 
nels. The weather conditions were various and resulting system 
temperatures were 350-500 K. 

The '^CO 3 -2 d ata for NGC1333 were published in 
iKnee & Sandelll ( |2000|) , and kindly made available by the au- 
thors. The spatial sampling of the NGC 1333 observations was 
10" rather than the 5" of the more recent observations. In order 
to match the other datasets as closely as possible for analysis, 
the maps were resampled onto a 5" grid. This was done by inter- 
polating with a 14" Gaussian truncated at 40". This results in a 
map with the same intensity (Tmb) scale but with 20" resolution. 

The sample (Table [TJ was selected before the analysis of the 
corresponding SCUBA data was complete, and therefore is bi- 
ased towards the obvious bright clusters and previously known 
protostars. In Fig. |2] we show the outflow sample (marked by 
stars) on a plot of luminosity vs. c ore mass for the entire pop- 
ulation of submm cores in Perseus ( iPaper Ij) . The sample is in- 
deed biased towards high luminosities with no outflow maps of 
sources below 0.25 Lq. The proportion of the total number of 
cores decreases with decreasing mass. However, the fraction of 
Lboi > O.25L0 cores studied remains fairly constant except for 
the very lowest mass sources below our completeness limit at 
10 K of Menv < 0.5 Mq (5cr). For the very highest mass sources, 
Menv > I2M0, all sources were observed as there were no low- 
luminosity sources in the SCUBA sample. Thus the most signif- 
icant bias of the sample is its lack of low luminosity sources. 

Small outflow maps provide a good estimate of the momen- 
tum flux (iFuller & Laddll2002l: iMangum et al.lll998h but, except 
for very small flows, do not provide estimates of the full extent, 
mass, kinetic energy, or dynamical timescale of outflows. Our 
map radius of 1' corresponds to 0.1 pc assuming a distance of 
320 pc for the Perseus molecular clouds, or a dynamical time of 
2000 years for an outflow advancing at 50 km/s. 



3. Outflow activity 

Maps of the outflows superimposed on the submm continuum 
map are shown in Fig. [1] We identified sources with outflows 
based on the objective criterion of '^CO linewings more than 
1.5 K (3cr) at 3 km s"' from the line centre (as measured in 
C'^O). Of the 51 sources in Table [1] observed in '^CO 3-2, 37 
sources (73%) were classified as outflows on this basis. Spectra 
for these sources at the source position and the outflow lobes are 
shown in Fig. |9] (online version only). In some cases the outflow 
lobe spectra are missing due to confusion or because the driv- 
ing source lies at the edge of the mapped area. From the maps 
(Fig.lTJ, we identify 13 cases where the broad linewings may be 
wholly or partly due to nearby flows driven by other sources in 
crowded regions. These uncertain flows are marked with 'y?' in 
Table[T] The remaining 14 sources show no evidence for outflow 
activity in the form of either line wings or red/blueshifted lobes, 
as shown in Fig. [TO] (online version only). 

All but four of the sources with outflows on the basis of 
linewings a lso have a near-infrared detection from Spitzer IRAC 
(iPaper iJ) . (|j0rgensen et al.l |2007|) . The exceptions are 17, 40, 
51 (SK16) and 53 (SK26), all of which are classified as star- 
less (Table [TJ. In all of these cases the linewings are likely due 
to confusion. In the following analysis of the outflow proper- 
ties, we leave out these four sources, assuming that they are in- 
deed starless. For most of the protostars it is the driving source 
which is being detected by Spitzer but in some cases the infrared 
emission may be due to HH objects (examples are HH21 1 (12), 
and (15)). In 9 of the 13 uncertain outflow sources, their proto- 
stellar nature is confirmed by IR detections. All but one of the 
sources with Spitzer IRAC detections is observed in CO have 
molecular outflows. The exception is in NGC1333: the Class 
source 61 lies to the north of the luminous Herbig Ae-Be star 
NGC1333 IRAS 8/SSV3 and B star IRAS 9/BD+30°549 and 
the interaction with radiation and winds from these sources may 
well explain the outflow non-detection. The detection of out- 
flows from all but one Spitzer-confirmed protostar shows that 
both types of observation are similarly eff'ective in identify- 
ing protostars among submm cores. Outflow-based identification 
avoids the Spitzer confusion with extragalactic objects but can 
be difficult in crowded regions where the outflows overlap. 

3. 1 . Outflow detection statistics 

If we consider the outflows as a function of core mass (Fig. [3] 
top), we find that the fraction of sources with outflows decreases 
with decreasing core mass. In order to estimate our ability to 
detect outflows in the lower mass bins, we performed the fol- 
lowing simulation. We considered the 14 outflow sources in the 
three highest mass bins (Menv > 1 Mq). We then simulated 
outflows from lower mass cores by scaling the outflow spec- 
tra (Fig. |9] online version only) for these 15 sources accord- 
ing to the modelled core mass as (intensity/original intensity) = 
(model mass/original mass)' The motivation for this power 
law scaling comes from Fig|4] which shows a rough correla- 
tion between momentum flux and core mass with a gradient of 
1.5 + 0.5 (see Sect. |4]i. We also formed the simulated outflow 
maps by multiplying the original (x, y, v) datacube by these scal- 
ing factors and adding back on the unsealed noise levels (from 
sampling the noise in the outer velocity channels of the unsealed 
maps). Several scaling factors were used, corresponding to mov- 
ing the core masses an integer number of mass bins to the left. By 
this means we generated an artificial population of spectra and 
maps of outflows for cores in the range O.6-5OM0. We tested 
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Fig. 1. i^co 3-2 outflow maps of IRAS 03292+3039 (76), IRAS 03282+3035 (77), B5 IRSl (78) and neighbouring submm source 
(79), IRAS03235+3004 (80), and IRAS 03271+3013 (81). Outflows are displayed as a renz ogram, equivalent to superimposing a 
series of channel maps each with one contour at a fixed intensity level (eg. iKregel et al.l2004h . Velocity ranges were as follows: (76) 
2.0 K km s-', -5-4 (blue) and 9-15 km s"' (red) ; (77) 3.0 K km s"', -5-5 (blue) and 7-15 km s"' (red); (78) and (79) 2.5 K km s"', 
-5-8 (blue) and 12-15 km s"' (red) ; (80) 2.0 K km s"', -5-3 (blue) and 6-15 km s"' (red); and (81) 1.0 K km s'', -5-6 (blue) and 
8-15 km s"' (red). Colourscale: 850yL/m SCUBA map with black contours at 100, 200, 400, 800, 1600 mJy/beam. The area observed 
in '^CO 3-2 is outlined in green. 
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B*i;B*3<J" lltj' VZ' Jfi* 



Fig. 1. continued. Region to the SE of IC348 The renzogram intensity level is 3 K km s ' and included channels from -5-7 km s ' 
(blue) and 1 1-15 km s ' (red). Region in the N of NGC1333. The renzogram intensity level is 2 K km s ' and included channels 
from -5-1 km s"' (blue) and 12-15 km s"' (red). 850/im contours and colourscale are as before. 
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Fig. 1. continued. IC348, region to the S of IC348, including HH211 (12) and IC348 mms (13). The renzogram intensity level is 
2 K km s ' and included channels from -5-6 km s ' (blue) and 1 1-15 km s ' (red). Bl ''CO 3-2 outflow maps of Bl, including 
Bl-c (1), Bl-bN/S (2), bl-d (4) and IRAS 03301+3057 (7). The renzogram intensity level is 3 K km s ' and included channels 
from -5-4 km s"' (blue) and 9-15 km s"' (red). 850jum contours and colourscale are as before. 
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Right Asoeftsicm (JSOOO) 



Fig.l. continued. L1448 including L1448 N/NW (27), L1448 C (28) and L1448 IRS2 (30). The renzogram intensity level is 
3.0 K km s ' and included channels from -5-0 km s ' (blue) and 8-15 km s ' (red). L1455 including RN015 (35), PP9 (37), and 
L1455 FIR2 (39). The renzogram intensity level is 1.8 K km s ' and included channels from -5-2 km s ' (blue) and 1 1-18 km s"' 
(red). 850yum contours and colourscale are as before. 
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Fig. 3. Top: Fraction of outflow sources as a function of core 
mass for sources with definite outflows (solid thick Une). The 
thin Une shows the number of Class I sources as classified in 
iPaper ll The ratio of protostars to starless cores decreases with 
decreasing mass. Bottom: As above with the fraction of proto- 
stars corrected for the outflow detection statistics in that mass 
bin. The greyscale gives the probability distribution of the frac- 
tion of protostars (see Appendix). 



the spectra according to our criterion for linewings (> 1.5 K at 
3 km s ' from line centre). We also additionally d isplayed these 
maps using our usual visualisation software kvis (iGoochll 19961) 
and considered whether we could visually detect the resulting 
outflows above the noise. In all cases where the linewings satis- 
fied the objective criteria, we could see the outflows in the maps. 
The detection statistics for our simulated population are given in 
Table |2] 

There is a decline in our detection ability at lower mass, and 
in the lowest mass bins we detect few of the model outflows. 
In the lowest two mass ranges in Fig. |3] 0.4-2.5 M©, and 1 
respectively out of 14 of the simulated low-mass outflows were 
detected. In the next highest mass bin (2.5-6 Mq), the detec- 
tion rate improved and 29% of simulated outflows were detected. 
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Table 2. Fraction of simulated outflows detected in each mass 
bin, compared to observed detection fraction. 



The mean detection rate below 6 Mgis predicted to be 12%, or 
roughly 1 in 8 outflows detected. 

The outflow detection rate predicted by the simulation is 
lower than the actual detection rate for all sources below 6 Mq. 
This can be seen in the lowest three mass bins of Table |2] This 
is true even if we include all the submm cores from the SCUBA 
survey, and not just those with existing outflow observations. We 
have already detected 14 outflows among 75 submm cores with 
masses below 6 Mq (Fig. |2]i, or 19% compared to an average 
of 12% for the simulated flows. This suggests that our assump- 
tions about how outflow linewings scale with core mass are pes- 
simistic, and we are more likely to detect flows than the simula- 
tion predicts - due to higher optical depth in more massive flows, 
for example. 

The bottom panel of Fig. [3] shows our best estimate of the 
true fraction of outflow sources in each mass bin, corrected for 
the detection statistics using (true outflow fraction) = (fraction 
of outflows detected) / (detection probability in this mass bin). 
The calculation of the uncertainties (greyscales) shown in the 
figure is given in Appendix lAl In mass bins where the observed 
detection rate equals or exceeds the simulated rate, this model 
would predict that all the apparently starless cores contain un- 
detected protostars. Only in the 6-16 Momass range is the pre- 
dicted number of outflow sources smaller than the number of 
apparently starless cores. If this model is correct, then we fail to 
detect outflows from the majority of protostars at masses < 6Mq. 

Thus, the modelling suggests that the low number of outflow 
detections from low mass cores is a selection effect. 

If this is the ca se then this outflow survey and Spitzer c2d 
lEvans et al. I (I2003h must have very similar sensitivities to em- 
bedded protostars as only one mid-IR detected sources in the 
sample had no corresponding outflow and conversely only four 
outflow sources had no mid-IR detection. Therefore any selec- 
tion effect must operate both in the infrared and molecular lines. 
Orientation could help produce such close agreement as infrared 
sources are hard to detect where disks are edge-on and outflows 
are hard to detect when they lie close to the plane of the sky. 
Assuming random outflow orientation, though, the number of 
infrared non-detections due to edge-on disks (within 10° of the 
line-of-sight) can only be ~ 20% of the total protostellar frac- 
tion. 

This modelling suggests that the low-mass cores without out- 
flow or Spitzer detections may also be protostellar, but of course 
does not prove that they are. For one thing, the assumptions 
in our model about how outflow masses vary with core mass 
may not be correct. Only further searching for outflows, infrared 
or radio sources towards the apparently starless cores can de- 
termine the true nature of the apparently starless cores. Note 
that such observations would need to be with higher sensitivity 
than this outflow s urvey, the Spitzer Cores to Disks programme 
(lEvans et al.l2003l) , or its successor, the Spitzer Gould's Belt sur- 
vey, in order to identify less luminous protostars. 



8 



Hatchell et al.: SCUBA Perseus survey - III. Outflows 



4. Momentum fluxes 

Outflow momentum fluxes (or outflow luminosities) are given 
in Table [T] for all the possible outflow sources. The method is 
developed from the momentum flux per beam calculation of 
[Fuller & Laddl (l2002h . rather than calculating the momentum 
flux within an annulus jBontemps et al.|[T996 ). though the two 
methods should yield comparable results. The momentum flux 
crossing a cylinder radius 0,. perpendicular to the plane of the 
sky and centred on the source, contributed by the mass within 
one beam, is calculated individually for each beam as: 



A^B(Vobs) Vpbs/ cos/ 

£>6'B/vobs tan/ 



where Vobs = \vi - vlsrI is the observed velocity offset from 
the core systemic velocity vlsr (determi ned from C'^O observa- 
tions of the same region at 1' resolution (iPaper Ih ) , / is the angle 
the outflow makes with the line-of-sight, D the source distance 
and 0B the beam FWHM. The numerator gives the momentum 
in the beam (v - Vobs/ cos / is the outflow radial velocity) and the 
denominator the time taken to cross the beamwidth at velocity v. 
The sum is over velocity channels. 

The mass within one beam Mb is calculated from the 
CO integrated intensity over the velocity chan nel Icq Av as- 
suming a CO abundance o f Xqo - 10 dPrerking et alj 
119821: I Wilson &Roodl 1 19941) and a LTE excitation at an 
outflow temperature of 50 K resulting in Nh-, - 2.5 x 
jq19(- Jt^^ dv/Kkm s"') cm"^. A temperature of 50 K is cho- 
sen as there is increasing evide nce that outflow s contain warm 
gasj Hatchell et al. 1999b a; Ni sini et al.l 120001: iGiannini et all 
1200 li) but the conversion is relatively insensitive to the tem- 
perature, increasing by only 40% over the temperature range 
20-100 K. The mass within one beam is calculated as Mb - 
(n/4)D^0lNn,mn,. 

The momentum flux contributions from individual beams are 
then summed in annuli, according to: 



2n9,. 



sm I 



where O,- is the angular distance of each pixel from the out- 
flow centre, D is the source distance, and the average is over 
pixels between radii 6i and 02 with A^pix the number of pixels be- 
tween those radii. As the contribution to the momentum flux is 
calculated individually for each beam, and varies little between 
neighbouring annuli, we set the inner radius to zero. The factor 
IndrlNOs corrects for the overcounting because the 14" beam 
size is larger than the 5" pixel size. 

The driving source positions were taken to be the SCUBA 
emission peaks from Paper I. To reduce the contribution due to 
emission outside the outflow region, the outflow areas to be in- 
cluded in the momentum flux calculation were selected by mask- 
ing with '^CO integrated intensity at a cutoff of 3 times the noise 
level in the integrated intensity map. We default to calculating 
the momentum flux within a circle of outer radius 30" unless it 
was clear from the map that the outflow was shorter than this or 
there was potential confusion within this radius, in which case 
the radius of the circle was reduced. In practice we find the 
choice of circle or annulus and the radius of the circle (except 



for the case of confusing flows) makes little difference to the re- 
sulting momentum fluxes, as expected at radii where there is no 
large variation in the injection of momentum. 

The inner velocity channels, contaminated by ambient emis- 
sion, were cut out to + 2 km s"' and the highest velocity consid- 
ered was ±15 km/s. Ideally, for a fair comparison of the momen- 
tum flux between strong (massive) and weak flows, one wants 
good signal-to-noise (S/N) in the linewings out to similar veloc- 
ities. With roughly constant sensitivity and line strengths vary- 
ing by an order of magnitude the strong linewings can be mea- 
sured with good S/N to higher velocities in the stronger flows. 
Adjusting the velocity range on a flow-by-flow basis according 
to where the S/N is good therefore introduces a strong bias in 
momentum flux towards massive flows. Not only are these more 
massive but also more of the high velocity channels are included 
in the calculation. Therefore we consider it fairer to use the same 
velocity cut for each source, even though this excludes the con- 
tribution from the highest velocity channels in all flows. 

The uncertainties on the momentum flux, listed in Table [1] 
are the root sum-of-squares of the uncertainties of each con- 
tributing pixel. The pixel uncertainty is taken to be the difference 
from the mean momentum flux in the contributing (integrated in- 
tensity > 3cr) region. 

We have not applied a c orrection factor for eit her inclina- 
tion or optical depth, unlike iBontemps et al.l ([1996'). Therefore 
our m omentum fluxes are strictly lower limits. Bontemps et ^ 
(1 19961) point out that inclination / reduces the momentum flux 
by a factor sin//cos^/ = 2.9 for a mean angle of 57.3 de- 
grees. For an optical depth t, uncorrected CO column densities 
are underestimated by a factor t/(1 - exp(-T)) or a typical fac- 
tor of 3.5 so that true momentum fluxes are typically a factor 
of 10 greater than measured. Outflow optical depths are high 
even in the J - li-l line with with T3_2 > 4 in the line wings 
(iMargulis & Ladalll985l:lHogerheiide et"ai1ll998l) . 

The (uncorrected) momentum flux results for all the possible 
outflow sources are given in Table [1] Comparing the sources in 
common, our momentum fluxes ar e typically lowe r by a factor 
of a few than previous analyses by iBontemps et al.l (Il996h (tak- 
ing the optical depth / inclination correction into account) and 
kneesandellOO. The main reason for this is that we assume the 
outflows are warm (50K), which redu ces mass estimates by a 
factor of ~ 4 over the lOK assumed bv IBontemps et all (119961) . 
This is a systematic difference, reducing the momentum fluxes 
for all sources by the same fraction. For flows with high velocity 
gas, we can also underestimate momentum fluxes because we do 
not include the highest velocity channels (as discussed above). 
For two sources, IRAS 03282+3035 (77) and B5 IRSl (78), we 
measure momentum fluxes more than a factor of 10 lower than 
previous estimates. We have no reason to suspect the calibra- 
tion for these two maps, and they do not have particularly broad 
linewings (as can be seen in Fig.|9] online). It is possible that in 
these flows, the momentum flux measured at greater radii in the 
flow is larger than in the region covered by our 2' maps. 

4.1. Momentum flux vs. Lboi, Menv and Tho\ 

The outflow momentum flux Fqo is plotted versus core mass 
Menv, bolometric luminosi ty Lboi and bolometric temperature 
Tboi in Fig. m (see iPaper III for the calculation of those quanti- 
ties from the SEDs). The well known rel ationship between mo- 
mentum flux and bolometric luminosity (IBontemps et al.|[T996l: 
Cabrit & Bertout 1992) is evident despite the intrinsic scatter 
in the momentum flux due to uncorrected outflow inclination. 
This correlation is often taken as evidence that the outflow and 
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Fig. 4. Relationships between momentum flux and (top to bot- 
tom) bolometric luminosity Lboi, core mass Menv, and bolometric 
temperature Tboi- Sources marked by open stars (blue) and aster- 
isks (red) are Class I and Class respectively based on the SED 
dPaper II ). The top two panels show best fits corresponding to 
Fco oc Lboi''^ and Fqo °^ A/env''^ respectively. Uncertainties on 
Fco (see Sect. |4]l are shown; uncertainties on Lboi are typically 
30%, Menv at least a factor 2, and Tboi a few Kelvin. 



the accretion luminosity are both being powered by the same 
mechanism, infall, and despite the details of the magnetic ejec- 
tion/accretion process in the inner few tens of AU of the disk, 
the energy input into each remains strongly linked. The correla- 
tion is tighter towards higher Lboi and holds for both Class and 
Class I sources (marked by red asterisks and blue open stars, 
res pectively; s ee Paper If for details of classification). 

iBontemps e t al. (1996J found a decline of momentum flux 
between the Class and Class I phases which has been taken 
as evidence that outflow power declines with age. In our sam- 
ple, this is less clear cut. The mean logio(/^co/A/okm s^'yr"') 
of the Class population is -6.0 with a sample standard devia- 
tion of 0.6, close to the Class I mean of -6.1 (standard deviation 
0.8). A histogram of the momentum flux distribution for Class 
and Class I sources (Fig.|5j shows that both of the distributions 
are very flat. Applying the Kolmogorov-Smirnov (K-S) test in- 
dicates a 27% probability that the two samples come from the 
same parent distribution and therefore there is no evidence for 
evolution in momentum flux between Class and Class I. There 
is considerable scatter with objects of both classes with both low 
and high-Fco- 

In particulai-, three Class I objects- NGC1333 SVS13 (43), 
NGC1333 ASR114 (45) and NGC1333 HH7-11 MMS4 (50)- 
have particularly high momentum flux and four Class sources 
- bl-c (1), bl-d (4), Bl SMM5 (76) and IC348 SMM3 (15) - 
have low momentum flux. The third panel of Fig|4] in which 
momentum flux is plotted against Tboi, shows clearly that the 
Class Is (on the right hand side of the plot) have more scatter 
in momentum flux than the Class Os. The two high Tboi sources 
are the well-known Class I sources L1455 PP9 (37), which may 
contain multiple protostars, and L1455 FIRl/2 (39). 

All three of the high momentum flux Class I cores, 
NGC1333 SSV13 (43), NGC1333 ASR114 (45) and HH7- 
11 MMS4 (50) are in the confused NGC1333 region. HH7- 
1 1 MMS4 (50) is a confused source where a major contribu- 
tion to the momentum flux comes from the HH7-11 driving 
source SSV13. Both SSV13 and ASRl 1 4 contain multiple pro- 
tostars as evidenced by radio emission dRodriguez et alj[l997l ; 
iKnee & S andell 2000). Can this explain their high momentum 
fluxes despite Class I classification? NGC1333 SSV13 (43) is 
a multiple source (SSV13 and 13B) potentially consisting of 
a Class I source (providing the infrared-based classification) 
and a Class source (powering the outflow). The flows from 




og^Q(Fj.Q/Wsun km s-1 yr-1) 



Fig. 5. Histograms of momentum flux for Class (solid) and 
Class I (dotted) protostars. 
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Fig. 6. Fco c/ Lboi (dimensionl ess) vs. Meny/Lb oi"'^ (in Mq nd 
respectively). iBontemps et aP (See also 1 19961 Fig. 7). Sources 
marked by open stars (blue) and asteri sks (red) are Class 1 and 
Class respectively based on the SED (iPaper ij) . 

SSV13/13B are confused and the momentum flux inc ludes con- 
tributions from both, though iKnee & Sandel (l2000l) attributed 
more momentum flux to the Class I SSV13 than the poten- 
tial Class counterpart SSV13B. Like wise ASR114 (45) con - 
tains multiple embedded radio sources dRodriguez et al ] ll997h . 
and the outflows are confused by the flow from the 4' dis- 
tant SSV13B which is believed to drive the nearby HH 12 
(iKnee & Sandelll 120001) . and the momentum fluxes could be 
overestimated by the contribution from the flow from the south. 
It is possible that in both these cases the high momentum flux is 
due to a Class source, either in the same submm core or nearby 
- an open question which could be resolved by interferometric 
observations of the flows. 

However, the Class sources with low momentum flux - bl- 
c (1), bl-d (4), IC348-SMM3 (15) and B1-SMM5 (76) - are less 
easy to explain away. Three of these sources (1,4, 76) also have 
low bolometric luminosity for their mass. Our sample, because 
of its submm-based sample selection, may well contain a pop- 
ulation of low-luminosity, low-Fco Class sources which were 
not included in earlier studies. 

There are physical reasons why outflow force may vary be- 
tween sources at the same evolutionary stage. The momentum 
in the driving wind depends on the mass accretion rate and in 
turn on the history of the mass reservoir available for accre- 
tion. For example, this may diff'er betwee n clustered and more 
isolated protostars i Henriksen et al.|[l997l) . though this does not 
appear to explain the results in Perseus. Early C lass source s 
may not have reached their full accretion rates ( ISmithll2000l) . 
The low momentum flux Class sources bl-c (1), bl-d (4) and 
B1-SMM5 (76) which also have low bolometric luminosity are 
likely to be early Class sources which have not yet reached 
their full accretion power. The fraction of the wind momentum 
which is transferred to the molecular outflow depends on entrain- 
ment of the ambient material, which may vary in efficiency. In 
the competitive accretion model, mass may continue to flow onto 
already-accreting envelopes and boost accretion at late stages, 
which would lead to a natural sc atter in accretion rates . Jet 
ejection is known to be episodic dRaga et al] 1 19901 12004|) but 
the swept-up molecular outflow smooths out any variability on 
tim escales of ~ 100 years, such as those which produce jet knots. 

iBontemps et alj ( 1 19961) show in their Fig. 4 a plot in which 
momentum flux and mass are compensated for source luminos- 



ity. In Fig.|6]we show the analogous figure for our dataset. For 
our sources the relationship between mass and luminosity is 
Lboi A^env''^^ dPaper III) so the abscissa becomes M^m / L^oi^'^ 
and the ordinate Fco c/Lboi- The luminosity correction separates 
the Class I sources from the Class Os, which tend to be more mas- 
sive for the same luminosity. This separation in mass is seen in 
Fig.|6] with the Class sources lying to the right of the diagram. 
It also takes out the effect of distance, though this is irrelevant 
for our Perseus-only sample. 

The luminosity correction has to be treated with caution be- 
cause it does not differentiate between luminosity and momen- 
tum flux evol ution. Evidence for luminosity evolution comes 
from iPaper ij, which shows that Class sources have, on av- 
erage, lower luminosities at a given mas s than Class I sources. 
Some evolutionary schemes (ISmithll2000l) interpret this as an in- 
crease in luminosity between the Class to Class I phases. On 
the luminosity-corrected diagram, decreasing luminosity can be 
mistaken for increasing momentum flux. 

Nonetheless, in the luminosity-corrected diagram, the 
FcQc/Lboi values for the Class sources are similar to their 
Class I counterparts, with no obvious evolution between the two 
stages. 

There is a general trend of increasing momentum flux with 
Menv (middle panel of Fig|4]i although there is more than an or- 
der of magnitude of scatter, as found for Fco-i'boi- A physical 
explanation for such scatter is that the energy in the mass infall, 
which powers both the mass outflow and the bolometric lumi- 
nosity, is dependent on the infall rate which is not set by t he cur- 
rent m ass of the core but by its past history. Hogerheii de et al.l 
(1 19981) found a better correlation of Fqo with Menv than Lboi, but 
we see no evidence for this. 

In summary, there is broad agreement between our results 
on the momentum flux and previous studies, with similar scat- 
ter on the correlations with Lboi or Menv despite the consistent 
derivation of these quantities and the single distance of our sam- 
ple. A difference is the significant population of low-Fco Class 
sources and the small number of high-Fco Class Is which sug- 
gest that evolution of Fco from Class to Class I is not a simple 
decline but that there may be variation within each class, eg. as 
luminosity and outflow power increase at the beginning of the 
Class phase. 



5. Mass loss rates and core destruction timescales 

A significant proportion of the momentum from the driving 
jet/wind is deposited in the protostellar core and removes matter 
from t he core in the molecular outflow (Lada 1985V Parker et al.l 
il99lt iFuller & Ladd 2002; Arce & Sargent 2006). The time it 
would take to clear a protostellar core of mass via the outflow is 
a constraint on the lifetime of the embedded protostellar phase, 
especially when taken together with an indication of whether the 
source is at the start (early Class 0) or end (Class I) of its embed- 
ded lifetime. 

In addition to the outflow, mass is lost from the protostellar 
core by accretion onto the disk/star. The ratio of mass ejected in 
the wind/jet to mass accretion, /, varies significantly with dif- 
ferent models for jet/wind la unching. Observatio ns suggest that 
/ =^ 0.3 is a good estimate (iRicheret al.ll2000l) . Although the 
mass loss in the jet is then less than half of the total accretion, 
momentum transfer to the slower molecular outflow can lead to 
mass loss rates in the outflow that are a factor of 10 higher than 
the total mass accretion rate: Mco = fM^ccivw/vco) - lOMacc, 
taking the jet ejection fraction /„ = 0.3 and typical values for the 
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jet and CO velocities of 150 and 5 km s ' respectively, and as- 
suming all the wind/jet momentum is deposited in the cloud. At 
early times when this is true, before the outflow has broken out 
of the cloud, then mass loss in the outflow dominates by an order 
of magnitude over the mass accreted onto the central object, and 
can be taken as a good estimate of the total mass loss rate. For 
more evolved flows less of the wind/jet momentum may be trans- 
ferred to the CO and in the most evolved sources the accretion 
rate may begin to dominate over the mass loss in the outflow. 

The outflow mass loss rate M can be calculated from the mo- 
mentum flux P by dividing by the mean outflow velocity v. v is 
a difficult quantity to estimate accurately. At low velocities it is 
difficult to separate outflow material from protostellar core. Any 
mass at high velocities can also make a large difference to v with 
correspondingly high uncertainties. Making an estimate using 
V = 2/ ^cokf-^LSRl/ YaI Icq where Iqo is the integrated spectrum 
over the observed inner section of the outflow, we calculate val- 
ues for individual flows which range from 1-10 km s ' depend- 
ing on where the integration limits - inner and outer velocity 
channel cutoffs - are set. We take a value of 5 km s ~' as typical, 
in agreement with previous estimates of P and M jZhang et alj 
120051: iKnee & Sandelll2000h . wifli an uncertainty of a factor of 
2. 

The resulting mass loss rates in the outflow range from 
1x1 0"^-6x 1 0"^ Mq yr" ' , from the momentum fluxes in Table[Tl 
and assuming a factor of 10 to take into account inc lination 
and op tical depth into account (see Sect. |4] and Bontem ps et alj 
( [19961) ). These core mass loss rates are similar to those cal- 
culated by [Fuller & Laddl ( |2002[) from the rarer CO isotopo- 
logues in molecular cores. Mass loss rates o f above 10'^ Mq yr 
have also been found for massive flows jZhang et al.1 120051: 
iBeuther et al. 2002; Churchwell 199^) and it is not a surprise 
to find rates approaching this in Perseus, which is forming stars 
up to spectral type B. 

The destruction timescale fa of the core at the current mass 
loss rate can be calculated by dividing the current core mass by 
the total mass loss rate, including the optical depth/inclination 
correction factor of 10. These timescales are plotted in Figs.|7] 
and range from 0.2-60 Myr, with a mean fd = 1 1 Myr These de- 
struction timescales are then lower limits to the lifetime remain- 
ing to each core, if outflows are the only source of mass loss and 
outflows only become weaker over time. This is believed to be 
true for all but the youngest Class sources. 

The sources with the longest lifetimes are the sources with 
the highest mass to momentum flux ratios. As luminosity and 
momentum flux are correlated (Fig.lUi, these are also the sources 
with high mass-to-luminosity ratios and cold temperatures, a 
signpost of an early Class source. Indeed, we find that the 
Class sources bl-c (1), bl-d (4), IRAS 03292+3039 (76) 
and the newly identified protostar IC348-SMM3 (15, Paper I, 
IWalawender et al. 2006) have long lifetimes (> 10 Myr). In all 
cases this is due to low momentum flux for the given core mass. 
Sources with particularly short destruction timescales ~ lO'* 
years are IRAS 03301-H3057 (7), Bl SMMll (10) and HH7- 
1 1 MMS (43). These sources have little mass in their cores and 
are believed to be Class I sources near the end of the embed- 
ded phase. Thus the objects with the most extreme lifetimes are 
Class for the longest lifetimes and Class I for the shortest life- 
times. The absence of short lifetimes for low sources is be- 
cause only high Menv sources can produce a low Tboi, and these 
tend to have a long lifetime. Low-mass sources have lower opti- 
cal depth in their cor es and hig her Tboi for the same evolutionary 
state, as discussed in lPaper 11 



The mean destruction timescale for Class Is (6 Myr) is 
shorter than for Class Os (14 Myr), suggesting that there might 
be some evolution between Class and Class I, with mass de- 
clining faster than outflow force. Yet the evidence for evolution 
is weak: the scatter in fa is large with both classes showing a 
range of destruction timescales from a few times 10^ to ~ 10^ 
years at the current outflow rate, as shown in the histogram of 
lifetimes (Fig[8l), and a K-S test confirms a 43% probability that 
these populations could be drawn from the same underlying dis- 
tribution. 

The mean lifetime of the protostellar stage is estimated to 



be 1.5-4 X 10^ years (Greene et al. 1994: 


Kenvon & Hartmann 


1995 


: IWilking et al.l 11989; .Kenvon et al.1 


1990; .Hatchefl et al. 



2006h . This is significantly shorter than our mean of 10^ years 
for the removal of the core mass by outflows. It was also recently 
suggested by Arce & Sargent (2006) that outflows may play an 
important role in mass loss and core evolution during the Class I 
phase. However, for the majority of sources. Class and Class I, 
either we have overestimated the timescales due to systematic 
errors or the bulk of the mass is removed by a means other than 
molecular outflows. This is compounded by the fact that not all 
of the mass accelerated in the low velocity flow will have ve- 
locities suflicient to remove it from the gravitational field of the 
central object, so the outflows then give an upper limit on the 
mass loss rate to the protostellar core. 

There are several factors which could introduce systematic 
errors in the calculation, which might result in an overestimate 
of the core destruction timescal es. We m ay have overestimated 
the core masses: as discussed in iPaper i5 our core masses lie at 
the high end of the range consistent with the data, and could be 
overestimated by a factor of 4 or more - though half of this is 
due to the uncertain distance, which affects dust and CO masses 
equally, so again it is hard to see how this can explain more than 
a factor of 2 difference between protostellar lifetimes and core 
destruction rates by outflows. The outflow mass loss rates are 
calculated using several uncertain factors in the conversion of 
CO 3-2 intensity to mass. Firstly, there is the optical depth and 
inclination correction, for which we assumed a factor of 10, and 
the outflow temperature, for which we assumed 50 K, though 
this cannot affect the masses by more than a factor of 1.5. The 
mean velocity is also uncertain to a factor 2 but the destruction 




0,5 

oq,Q(lifetime / Myr) 



Fig. 8. Histograms of core destruction timescales for Class 
(solid) and Class I (dotted) protostars. Sources with upper limits 
on the lifetime are excluded. 
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timescales would only reduce if it were lower than the assumed 
5 km s"'. 

If protostellar cores are not destroyed by outflows and ac- 
cretion, where does the mass go? A possibility is that we are 
detecting material that will ultimately form not single stars but 
small clusters, and stars will continue to form until t he mass is 
used up, as is seen in hydrodynamical models (eg. iBate et all 
l2003b . Many of the most massive cores ar e already known to 
contain multiple protostars (Sect l4lPaper Ilh . A further possibil- 
ity is that some cores will ultimately form massive stars and the 
molecular cores will be destroyed by ionisation. Finally, as mass 
is lost from the system and the core heats up some material may 
become gravitationally unbound and dissipate. 



6. Summary and conclusions 

We have searched for molecular outflows towards 5 1 of the 103 
known subm m cores in Perseus list ed in Paper III and based on 
the SCUBA dPaper J) and Bolocam jEnoch et al.ll2006l) surveys. 
The outflow sample is biased towards the higher luminosity 
submm cores but the fraction of cores sampled is fairly inde- 
pendent of mass except at the lowest masses. The main results 
of this study are the foUowing: 



1 . We detected broad linewings from 37/5 1 sources surveyed. 
Of these, 13 cases could be due to confusion with neighbour- 
ing flows, though 9 are confirmed protostellar by Spitzer in- 
frared detections. Only one source with a Spitzer detection 
has no apparent outflow. 

2. The well-known correlation of momentum flux with bolo- 
metric luminosity is present, though with considerable scat- 
ter. Momentum fluxes are not significantly higher for Class 
sources than Class I: on the basis of a K-S test, we cannot 
reject that the two samples are drawn from the same distri- 
bution. There are three Class I sources with apparently high 
momentum fluxes, possibly due to confusion with Class 
flows, and there are several Class sources with low mo- 
mentum fluxes, possibly because they are very young and 
have not yet reached full power 

3. Core destruction timescales based on constant outflow mass 
loss rates and core masses from submm emission range from 
10^-10** years, with no clear evolution between Class and 
Class I sources. These timescales are longer than the esti- 
mated lifetimes for protostars of 1.5-4 x 10^ years. There 
are several possible explanations for this discrepancy. These 
include the possibility that outflows may not be the primary 
mechanism which destroys the cores. Alternatively, although 
it appears unlikely that the core mass is overestimated due to 
systematic issues by a large enough factor to account for the 
timescale discrepancy, it is possible that we are overestimat- 
ing the fraction of the core mass bound to each protostar and 
that the material fragments forming multiple sources. 
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Fig. 1. continued. NGC1333 from data first published in iKnee &Sandel]l (120001) . The renzogram intensity level is 3 K km s ' and 
included channels from -5-3 km s ' (blue) and 12-18 km s"' (red). 850jum contours and colourscale are as before. 
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Table 1. Sources with outflow observations and resulting momentum fluxes. For each source we list the core number and name (see 
iHatchell et al.l ( l2006l) for positions and core properties), the outflow status; the core ambient velocity vlsr, derived from C"^0, the 
maximum radius used in the momentum flux calculation, 02 {Oi was set to zero); and the resulting momentum flux and uncertainty. 
Of the 51 outflow sources, 24 have definite outflows, 14 have no outflow detection, and 13 (9 with infrared detections) have possible 
but confused outflows. 



Nr. 


Name' 


Outflow 

(y/n?) 


Class- 


I'LSR 

km s"' 


02 


Fco 

10"' Mokm s"' yr"' 


1 


bl-c 


y 





6.4 


12.0 


1.7(0.8) 


2 


bl-bS 


y 





6.8 


30.0 


12.3 ( 1.0) 


4 


bl-d 


y 





6.7 


15.0 


1.8(0.3) 


7 


IRAS 03301+3057 


y 


I 


6.6 


30.0 


12.6(0.8) 


10 


Bl SMMll 


y 


1 


6.6 


30.0 


4.4 ( 0.4) 


12 


HH211 


y 





8.8 


30.0 


6.2 ( 0.8) 


13 


1C348 MMS 


y 





8.6 


30.0 


5.1 (0.5) 


14 




y? 


1 


9.1 


10.0 


1.2(0.3) 


15 


1C348 SMM3 
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' Common name. For positions, further na mes, refer ences and further data see |Paper"l3 . 
^ From SED-based evolutionary indicators jPaper I]) . 
^ Sources of potential confusion. 
At edge of mapped area. Possible S blue lobe. 
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Fig. 2. Lboi vs. Menv for all submm cores in Perseus dPaper it [Paper ij) with '^CO targets marked by small black circles and outflow 
detections additionally marked by large black circles. Class sources are marked as asterisks (*, red) and Class I by open stars 
(blue). The remaining sources (dots) are assumed starless. The horizontal line at 0.25 Lq marks the lower luminosity limit for the 
outflow sample, and the vertical dashed line at 0.3 Mq marks the 3cr detection limit for the SCUBA survey. 
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Fig. 7. C ore destruction times cales at current outflow mass loss rate, vs. left bolometric temperature; right current core mass. Class I 
sources jHatchell et al.ll20d6l based on the SED), are marked by (blue) open stars and Class sources by (red) asterisks. 
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Appendix A: Calculation of outflow fraction 
probability 

Two measured quantities ai'e involved in calculating the fraction of protostars 
(outflow sources) in each mass bin. The first is the observed fraction of outflow 
sources in each mass bin m. B = bjn (n is the total observed sources in the bin). 
The second is the detection probability in that inass bin calculated froin the sim- 
ulation: A = a/N (N is the total simulated sources in the bin). Both a and b are 
expected to follow a binomial distribution so we can easily calculate the corre- 
sponding uncertainties on B and C: they are just Bin (a, n. A) and Bin(fc, n,B) 
respectively, where Bin (k, n, p) is the binomial distribution for k successes from 
n trials with probability p. With our estimates for the detected outflow fraction 
A and detection probability B we can then estimate the true protostellar fraction 
in the bin C = B/A. 

Given that a and b are binomially distributed, we can also calculate the 
probability distribution for C = A/B by integrating the probability space P(A, B) 
where B = A/C is satisfied. 



'''' - // 



P(A,B)S(C -A/B)dAdB 




P(A) P(B = A/C)dA 



Bin (a, n, A)ln X Bin (b, n, B)IN dA 



This we integrated numerically for a range of values of C for each mass bin 
and the resulting probability distributions are plotted as greyscale in Fig. [3] 
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Fig. 9. '^CO 3-2 spectra for outflow sources showing linewings at the 1.5 K (3cr) level at least 3 km s"' from the ambient velocity 
(marked by inner solid and vertical lines). The outer vertical lines mark the range of velocities included in the momentum flux 
calculation. 




Fig. 9. continued. 
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Fig. 9. continued. 
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9. continued. 
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9. continued. 




Fig. 9. continued. 
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Fig. 10. '^CO 3-2 spectra for sources with no outflow, showing no linewings above the 1.5 K (3cr) level (upper horizontal solid 
line) at 3 km s ' from the ambient velocity (inner solid and vertical lines). 



